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ABSTRACT 

The aim of the present paper is to propose an algorithm 
for a new ODE-solver which should improve the abilities of 
current solvers to handle second order differential equations. 
The paper provides also a theoretical result revealing the 
relationship between the change of coordinates, that maps 
the generic equation to a given target equation, and the 
symmetry 7?-groupoid of this target. 

Categories and Subject Descriptors 

I.1.2H.4 [Symbolic and Algebraic Manipulation]: Com- 
puting Methodologies — Algorithms 

General Terms 

Algorithms 

Keywords 
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1. INTRODUCTION 

Current ODE-solvers make use of a combination of sym- 
metry methods and classification methods. Classification 
methods are used when the ODE matches a recognizable 
pattern (that is, for which a solving method is already imple- 
mented), and symmetry methods are reserved for the non- 
classifiable cases - Fig. [1] Using symmetry methods, the 
solvers first look for the generators of 1-parameter symme- 
try groups of the given ODE, and then use this information 
to integrate it, or at least reduce its order 01 [5] • 

In practice, present solvers are often unable to return 
closed form solution. Consider for instance the following 
equation 
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Figure 1: General flowchart of typical ODE-solver. 



which admits only one f -parameter symmetry group. Using 
this information, actual solvers return a complicated first 
order ODE and a quadrature. Clearly, such output is quite 
useless for practical applications. More dramatically, con- 
sider the following equation 



y = 



(2) 



When applied actual solvers, output no result. This failure is 
due to the fact that the above equation does not match any 
recognizable pattern and has zerodimensional point symme- 
try group(oid). Thus neither symmetry method nor classi- 
fication method works. 

Our solver (the implementation is in progress) is designed 
to handle such equations. It returns an equation from Kamke's 
book [S], equivalent to the equation to be solved, and the 
equivalence transformation ip. Thus, for the equation ((T| we 
obtain the Rayleigh equation y" -\- j/'* -\- y — Q (number 72 
in [S]) and the change of coordinates : {x,y) — > (x, y^/2). 
For the equation ([2]), we obtain the first Painleve equation 
y" = 6y^ + X (number 3 in igj ) and the change of coordinates 

: {x,y) (1/a;, y)- It is worth noticing that this transfor- 



mation can be composed with the point symmetries of the 
Painleve equation given by {x,y) (A^a;, Ay) with A^ = 1. 

To summarize the theoretical result of the paper, let £f 
denote the generic ODE y" = f(x, y, y') and $ an allowed 
O-groupoid acting on the variables x and y. Suppose that 
<I> is given by quasi-linear Lie defining equations. Define 
Sf := aut(£/) n"I> where aut(£/) denotes the (local) contact 
symmetry 23-groupoid of £f. Given / £ Q(a;,j/,2/') and as- 
sume that the symmetry ©-groupoid Sf is zerodimensional. 
The transformation if £ $ mapping the generic equation £f 
to the target equation f / is called the necessary form of the 
change of coordinates. We shall see that this transformation 
exists and belongs to the differential field Q(/), endowed 
with the partial derivations {dx, dy,dyi), for almost any an- 
alytic function / satisfying (y3«(f/) = £j. In other words, 
1/5 is an algebraic function in / and its partial derivatives 
and can be obtained without solving differential equations. 
Moreover, the degree of ip is equal to card(iSj) which is finite. 
Note that, the use of the O-groupoids formalism is dictated 
by the non global invertibility of the transformation ip. 

As we shall see, ip can be obtained using differential elimi- 
nation. Unfortunately, such approach is rarely effective due 
to expressions swell. For this reason, we propose in sec- 
tion |4] a new method to precompute the transformation if in 
terms of differential invariants, for each target equation £j 
in Kamke's list. These invariants are provided by Cartan's 
method. In the last section, we present our solver. This 
solver uses 7 possible types of transformations $i, • • • , $7. 
Using Lie infinitesimal method we precalculate to each tar- 
get equation a signature. That is, the dimensions of the 7 
symmetry groupoids associated to the 7 groupoids $1 , ■ ■ ■ ,^7. 
If two differential equations are equivalent then their signa- 
tures match. This fact significantly restricts the space of 
research in kamke's list at the run-time (when the input 
equation £f is known). 

2. EQUIVALENCE PROBLEMS 
2.1 Groupoids 

Definition 1 (Groupoid). A groupoid is a category 
in which every arrow is invertible. 

Let {Q, o, s, t) be a category. Each arrow if G Q admits 
a source s(i/j) and a target t{ip) which are objects of this 
category. For each arrows a, P £ Q such that s(/3) = t(a), 
there exists a unique arrow P o a £ Q with the source s(a) 
and the target t(/3). If 5 is a groupoid, for each arrow if £ G, 
there exists a unique inverse arrow ip~^ such that ip~-^ oif = 
Ids(^) and ip o if~^ = Idt(^). 

Let X and U be two manifolds and x £ X. The Taylor 
series up to order q (i.e. the jet of order q) of a function 
f : X U , oi class C, is denoted j%f. The Taylor series 
of / about x is denoted j^f or f. We shall say that x € X 
is the source and /(x) £U \s the target of the g-jet jj/. 

Example 1. For instance, when X — U — C, we have 



jlf:= [x,f{x),f'{x),...J^'Hx)) 



-9-1-2 



For each q £N and each x G X, we set 3i{X, U) := IJ/ Af 
and ^(X.f/) := Uxgx Jx(^, f/)- We denote by J^(X,X) 
the submanifold of J'(X, X) formed by the invertible jets. 
J»(X, X) is a groupoid for the composition of Taylor 
series up to order q according to 



A{9°f) = {ff^^)9)ohl.f)^ 



(3) 



By definition, a ©-groupoid tl2j Q C J^(X, X) is a sub- 
groupoid of J^(X, X) formed by the Taylor series solutions 
(see def. |8]) of an algebraic PDE system called the Lie defin- 
ing equations. This system contains an inequation which ex- 
presses the invertibility of the jets. The set of smooth func- 
tions ip : X ^ X that are local solutions of the Lie defining 
equations of 5 is a pseudo-group denoted by VQ. We define 
dimC/ := dimC and, if dimC = 0, degC/ := degC where C 
is a characteristic set (see sect. |3]) of the Lie defining equa- 
tions. We have degQ = card(r0). 

Example 2 ($3). Let$3 (see table 1) be theV-groupoid 
of infinite jets of transformations (x,y) = ip{x,y) where 

X = X -\- C and y = ri{x, y). (4) 

The constant C £ C and the function 77 : — + C are 
arbitrary. <I>3 C J^(C^,C'^) is an infinite dimensional V- 
groupoid where the corresponding Lie defining equations are 



Xx = 1, Xy = 0, yy / 0. 



(5) 



Definition 2 (CJ-Invariant). An invariant of the V- 
groupoid Q C JJ°(X, X) is a function / : X — > C which is 
constant on the orbits of Q . 

Clearly, the sum, the product and the ratio of two invari- 
ant functions is still an invariant function. Consequently, 
invariant functions of Q define a field. 

2.2 Differential equations and diffieties 

Let £f denotes the generic ODE 



y 



= f{x, y, y , 



(6) 

Let M := J"(C, C) be the n-th order jets space [16] of func- 
tions from C to C. Let x := (a:, y, yi, . . . , y-n) £ C" be a 
local coordinates system over M where m := n+2 = dim M. 

Every differential equation £f defines a diffiety [11]. This 
diffiety is given by the manifold M and a set of 1-forms, 
called contact forms, satisfying the Frobenius condition of 
complete integrability. Contact forms are linear combina- 
tions of the basic contact 1-forms dy~yidx, dyi—y2dx, • ■ • , 
dyn — /(x)da;. Vector fields which are orthogonal to the con- 
tact forms are colinear to the Cartan field 

_ 9 d d CI \ d 

ox dy dyi dyn 

They generate a distribution denoted by A/ . A local isomor- 
phism If between two diffieties £f and £f is, by definition, a 
local diffeomorphism if : M ^ M such that 



</^.(A^ 



This jet is said to be invertible if f'{x) 7^ 0. The jet of the 
function Id about the point x is (x, a;, 1, 0, . . . , 0). 



2.3 Equivalence problem and symmetries 

Definitions (Equivalence problem), inequivalence 
problem (EPB) ts an ordered pair (M, $) where M = J"(C, C) 
and $ C JJ°(C^, C^) is a V-groupoid of point transforma- 
tions from to . 



There exists a unique prolongation of denoted <1>^"', that 
acts on M (see section 13.4. ip . Two differential equations 
£f and £j are said to be equivalent if there exists a local 
transformation tp : M ^ M satisfying the differential system 



Aj = 'y3*(A/) and <^ £ r$ 



(n) 



(7) 



The second condition means that ip fulfills the Lie defining 
equations of the 2?-groupoid 'l?'-"^ 

The system (O is fundamental and we shall see that it can 
be treated by two different approaches : brute-force method 
based on differential algebra (section [3)l and geometric ap- 
proach relying on Cartan's theory of exterior differential sys- 
tems (section U). It is classically known that the existence of 
at least one transformation tp can be checked by computing 
the integrability conditions of the system (0 , which is com- 
pletely algorithmic whenever / and / are explicitly known 
[51 1161 fT|. However, there is no general algorithm for com- 
puting closed form of ip. In the sequel, we shall show that if 
the function / is fixed such that a certain D-groupoid <Sj is 
zerodimensional, then Lp is obtained without integrating any 
differential equation. 

Definition 4 (Sf). To any EPB, with fixed target equa- 
tions f, we associate the'D-groupoid Sf C {M, M) formed 
by the Taylor senes solutions of the self-equivalence problem 



Aj = (T*(Aj) and a G 



(n) 



(8) 



Example 3. Consider the EPB (J^(C,C),$3) and the 
Emden-Fowler equation Ejr (number 11 in f^) 



y = — 2- 



(9) 



The Lie defining equations of the T>-groupoid Sj are given 
by the characteristic set 



- yp -3 3 - 

p = — , y =y , x = x 

y 



(10) 



This PDE system is particular. Indeed, it contains only non 
differential equations. We have AimSj — and deg<Sj — 3. 
We deduce that its associated pseudo-group 

r<Sj= {{x,y,p) -> {x,Xy,\p) I A"" = 1} 

is, actually, a group with 3 elements. 

Equivalence problem and associated V-groupoid 

Let X := J°°(M, C). Any EPB (A/, $) defines a C-groupoid 
g C ■JT{X,X) formed by the set of triplets 

(jx/, JxV, JvmI) 

where x € Af and the functions (/, ip, f) are local solu- 
tions of the differential system (0. The source of a triplet 
is the infinite jet j^f £ X and the target is the infinite jet 
J<pM I The composition of two triplets (j^/, j^ipx, j^^ fi ) 
and (jxi/i, J^if2, JX2/2) is the triplet {j^f, j^ip, 7x2/2) where 
we have ip := ip2 o ipi. 

Definition 5 (Specialized invariant). For each Q- 
invariant I and each function f : M ^ C, we define the 
specialized invariant I\f] : M —* £, by 



3. USING DIFFERENTIAL ALGEBRA 

The aim of this section is to use differential elimination 
to solve the EPB when the target function / is a Q-rational 
function, explicitly known and the I>-groupoid of symme- 
tries Sj is zerodimensional. 

3.1 The vocabulary 

The reader is assumed to be familiar with the basic no- 
tions and notations of differential algebra. Reference books 
are ^ and [10]. We also refer to [1 [8| [l]. Let f/ = 
{ui, ■ • • ,Un\ be a set of differential indeterminates. k is a 
differential field of characteristic zero endowed with the set 
of derivations A = ■ ■ • , dp\. The monoid of derivations 



9:= {9^9° 







(12) 



/[/](x) := /(j,/), xGA^. 



(11) 



acts freely on the alphabet U and defines a new (infinite) 
alphabet QU . The differential ring of the polynomials built 
over QU with coefficients in k is denoted J? = k{C/}. Fix an 
admissible ranking over QU . For f £ R, ld(/) € QU denotes 
the leader (main variable). If £ R denotes the initial of / 
and Sf £ R denotes the separant of /. Recall that Sf — ^ 
where v = ld(/). Let C C -R be a finite set of differential 
polynomials. Denote by [C] the differential ideal generated 
by C and by /[Cj the radical of [C]. Let He ■= {If \ 
f £ C}U {Sf I / £ C}. As usual, fulLrem is the Ritt fuU 
reduction algorithm JjOJ. If r =fulLrem {f,C) then 3h £ 
He : hf — r mod [C] . Then the normal form is defined by 
normaLform(/) := r/h. 

Definition 6 (Characteristic set). The set C C R 
is said to be a characteristic set of the differential ideal 



(1) C is autoreduced, 

(2) f £ c if and only i/ fulLrem (/, C) = 0. 

Definition 7 (Quasi-linear characteristic set). 
The characteristic set C £ R is said to he quasi-linear if for 
each f £ C we have deg(/, w) = 1 where v is the leader of f. 

Proposition 1. When the characteristic set C is quasi- 
linear, the differential ideal c := : C R is prime. 

3.2 Taylor series solutions space 

Let k :— Q(a::i, ■ ■ ■ ,Xp) be the differential field of coeffi- 
cients endowed with the set of derivations < -jS- , ■ ■ ■ , tt- f ■ 

\^ axi ' ' dxp J 

Let C be a characteristic set of a prime differential ideal 
c £ R. We associate to C the system 

(C = 0, ffc / 0) (13) 

of equations f — 0, f £ C and inequations h ^ 0, h £ He. 

Definitions (Taylor series solution). ^4 Taylor se- 
ries solution of the PDE system 1113}) above is a morphism 
/i : -R — > C of (non differential) Q-algebras such that 

[C] C ker and He n ker ^ = 0. 

The morphism /x defines an infinite jet where the source is 
s(/x) := {fj,{xi), . . . , fi{xp)) £ and the target is t(/i) := 
{fj,{ui), . . . , ^{un)) £ C". Thus, a Taylor series is a C-point 
of an algebraic quasi-affine variety. Its Zarisky cloture is 
an affine variety defined by the ideal c. The dimension of 
the solutions space of (|13|l is the number of arbitrary con- 
stants appearing in the Taylor series solutions fi when the 



source point x := s(^) £ C is determined. Let K be the 
fractions field Frac(7?/c). Recall that the transcendence de- 
gree of a field extension K/'k is the greatest number of el- 
ements in K which are k-algebraically independent. The 
degree [K : k] is the dimension of as a k- vector space. 
When tr deg(A"/k) = 0, the field K is algebraic over k and 
[-ftT : k] < oo. If f £ C, we denote rank(/) := {v,d) where 
V := Id / and d := deg(/, v). Let 

rankC — {rank(/) | / G C} 

IdC - {ld(/)|/GC} 
dimC — card(e[/\e(ldC)) 

degC H'i^Said/). 

fee 

Proposition 2. dimC = tr deg(_K'/k) ts the dimension 
of the solutions space of f //dimC — then the cardinal 
of the solutions space is finite and equal to deg C = [-K" : k] . 

3.3 Differential elimination 

Let U = Ui U U2 he a, partition of the alphabet U. A 
ranking which eliminates the indeterminates of U2 is such 
that 

Vi;i G et/i, Vu2 G et/a, V2yvi. (14) 

Assume that C is a characteristic set of the prime differ- 
ential ideal c = \/[C] : Hq' w.r.t. the elimination ranking 
QU2 >- QUi. Let Ri ~ k{Ui} be the differential polyno- 
mials fc-algebra generated by the set Ui. Consider the set 
Ci :— C n Ri and the differential ideal ci := c n -Ri. 

Proposition 3. Ci is a characteristic set o/ci. 

Consider the differential field of fractions K ~ Frac(i?/c) 
and denote hy a : R ^ K the canonical fc-algebra morphism. 
Let Ki be the differential subfield of K generated by the 
set a(-Ri). Then Ki is the fraction field associated to the 
prime differential ideal ci := c n The partition of the 
characteristic set 

C = Ci U Ca (i.e. C2 := C \ Ci). (15) 

enables us to study the field extension K/Ki. 

Proposition 4. tr deg{K/Ki) = dimCa. //dimCa = 
then [K : Ki]= degCa. 

3.4 The system © revisited 

3.4.1 Prolongation algorithm 

Our aim, here, is to prolong the action of "l? C JJ°(C^, C^) 
on the manifold M := J"(C,C). For each integer g > 0, 
define 

k(«) := Q{x,y,yi,...,y,) 
R^"^ := k(''{S,y,yi,...,y,} 

The differential field k'^''' is the coefficients field of the ring 
of differential polynomials -R*-'' endowed with the set of 

derivations | ^ , ^ , . . . , -^j- 1 . Let us assume that the Lie 
defining equations of "1> are given by a characteristic set 
(7(0) ^ ^(0) rpj^g D-groupoid acting on J-^ and pro- 
longing the action of $ is characterized by a characteristic 



set C^'-* C -R^'-* . The prolongation formulae [15] of the point 
transformation {x,y) {^{x,y), r]{x,y)) are of the form 

where y — ri{x, y) if q — 0. The computation of the char- 
acteristic set C*-'' is done incrementally using the infinite 
Cartan field ~ £+yi-§^+ y2-^ + ■■■ 

rig := D^riq-i ■ {D^O^^ 
(jii) ■- c^i~i) u |y, - normalJorm (y^, C^'-^^) | 

Proposition 5. // C^°^ is a quasi-linear characteristic 
set of $ then C^"'^ is a quasi-linear characteristic set of 
^(q) yj^j-^t^ tfig elimination ranking Qyq >- y ■ ■ ■ y 

e{y,x}. 

The previous proposition gives an efficient method to pro- 
long a O-groupoid $ without the explicit knowledge of trans- 
formations. 

3.4.2 EPB with fixed target 

Let us compute a characteristic set C[f] C R^"^{f} for the 
PDE system ^ where / is fixed Q-rational function. First, 
prolong C'"' up to the order n -I- 1 as above. Then C[f] is 
obtained by substituting in c'""'"^' the indeterminate yn+i 
by the symbol / and the indeterminate i/n+i by the rational 
function f{x, ■ ■ ■ , y„). 

Example 4. For the EPB (J^(C,C),$3), we have 
P-Vx- PVy = 0, 
yxx + 2py^y +p^yyy + f Vy - fix, y, p) = 0, (16) 

- 1 = 0, = 0, Xp = 0, J7p = 0, ?/y / 0. 

These equations constitute a quasi-linear characteristic set 
w.r.t. the elimination ranking Qf >- Qp >- Oy >- Qx. 
Hence, the associated differential ideal is prime. 

Corollary 1. The PDE system ^ (where f is a fixed 
Q-rational function) is a quasi-linear characteristic setC[f] C 
/?'"•'{/} w.r.t the elimination ranking Of >- Qyn >■■■>- 
e{y,x}. 

3.5 Brute-force method 

Using Rosenfeld-Grobner we compute a new charac- 
teristic set C[f] of the PDE system ([TJ w.r.t. the new rank- 
ing 0{j/n, ■ ■ ■ , y\,x} y 6{/}. We make the partition of 
C ~ C[f] as in ^ 

C = CfU (17) 

where C/ := Cnk("){/} and C^~C\Cf. 

Proposition 6. The transformation tp does exist for al- 
most any function f satisfying the PDE system associated 
to the characteristic set C/[/]. The function x = fix) is 
solution of the PDE system associated to C^[f]. 

If dim [f] = Q, one can calculate ip by an algebraic process 
without integrating differential equations. 

Definition 9. When dimC,^[/] = 0, the algebraic sys- 
tem associated to C,p[f] is called the necessary form of the 
change of coordinates x = v'(^) • 



Example 5. Consider the EPB {J^{C,C), ^3). Suppose 
that the target £j is the Airy equation 

_// 

y = xy. 

In this case, Rosenfeld-Grobner returns C^[f] and C/[/] 
resp. given by il8d and U9d 



y^x = -fyy + pfpVy - 1/2 p'^fppyy + yfy - 1/2 y fa: 

-1/2 yfppf + 1/4 yfp^ - 1/2 ypfyp 
Vxy = -1/2 fpyy + l/2pfp^pyy 
Vyy ~ ~l/2/ppJ/i/, Vp = 0, 
X = /, -l/2/,j,-l/2/pp/ + l/4//-l/2p/^p 



(18) 



fxxp = 2 fxy + fpfxp — 2 + p' fyyp ^ fppfx + 

fxyp = 2 fyy — pfyy' 

fxpp = fyp ~ Pfypp 

fppp — . 



(19) 



We have dim [/] = 3 which means that the transformation 
X = i^(x), when f satisfies Cf[f], depends on 3 arbitrary 
constants. 



Theorem 1. The following conditions are equivalent 

(1) dim(C^[/]) = 0, 

(2) dim(5j) = 0, 

(3) deg(5j) < 00. 

In this case, degiSj = deg(C,p[/]) = degiys. 

Proof. Define 

Gf := {(jx/, JxV', Jx/) G e I / determined} . 

5/ is an algebraic covering of M defined by the characteristic 
set C^[/]. The X>-groupoid Sj C Qf is defined by diff'erential 
system (|8]) i.e. the characteristic set Cct[/]. Figure |2] shows 
that Sf acts simply transitively on Qj. 




3xq f 



Jx/ 



Figure 2: Simply transitive action of Sf on Qj where xq 
<^o(x) and x = ip{x) 



Example 6. Consider the EPB (J^(C,C),$i) where <l>i 
is defined in table 1. Assume that the target equation £j is 



-II -3 

y =y ■ 



Here, Rosenfeld-Grobner returns dplf] and Cf[f] resp. 
given by (2^ and l[2l\) 



(20) 



f = 1/12 (i fy - 2 fxp- 2 fppf-2pfyp + fp' 



iA fy fxp 2 fppf 2pfyp-\-fp ) X 



fxxyp (4 fy 2 fxp 2 fpp.f 2pfyp fp ) X 



(4 /y — 2 fxp ~ 2 fppf — 2pfyp + fp"^) ^ X 



(21) 



fxpp — fyp pfypp 
fppp ~ 0. 

Consequently dimC,p[/] = and degC,p[/] = 2. Thus, ip 
is the algebraic transformation of degree 2, given by equa- 
tions ESI). 



3.6 Discrete symmetries o-groupoids 

The self-equivalence problem, is in fact, the EPB when 
the PDE system ((7} is specialized by substituting the symbol 
/ by the value /(x), that is 



/ := f{x, y, 



(22) 



After specialization, the difi'erential system Cf[f] constrain- 
ing the function / is automatically satisfied (since there ex- 
ists at least one solution x = (t(x) of the problem, namely 
a = Id). The symmetries a are solutions of a characteristic 
set Ccr[.f] obtained form C^lf] by the specialization (|22p . 

By definition, the degree of an algebraic transformation 
X — fix) is the generic number of points x when x is deter- 
mined. 



Choose a point xo in AI. For every <^o G r$, define the 
rational transformation Sj Qj 

Jro^ ~* (Jxo'^) ° (JxV^o), (o- G rSf). 

In fact, according to the Taylor series composition formulae, 
this transformation is birational. Thus, the one-to-one cor- 
respondence between the two algebraic varieties Qf and Sf 
is birational. Consequently, the two characteristic sets [/] 
and Cct [/] have the same dimension and the same degree. □ 

Lemma 1. The rank of the characteristic set Cip[f] is sta- 
ble under the specialization II22V i. e. rank [/] = rank Co- [/] • 



Proof. The specialization (|22|) transforms the character- 
istic set Ciplf] to Ccr[f]. A fall of the rank of C,p[f] during 
the specialization contradicts the existence of birational cor- 
respondence between Qjr and Sf. □ 

Remark 1. When the transformation x := <p{x) is locally 
bijective but not globally, Sf and Sf need not to have the 
same degree. Indeed, consider again the groupoid $3 and 
the equations 



y 



&y^ + x~ 2y'' 
2^ 



and y" — 6y -'r x 



which are equivalent under {x = x, y = y^). The corre- 
sponding symmetry group are respectively given by 

TSf = {{x, y) ^ {x, \y) | = 1} and VSf = {Id}. 

They have the same dimension but different cardinal. 

3.7 Expression swell 

In practice, the above brute-force method, which consists 
of applying Rosenfeld-Grobner to the PDE system ((Tjl, 
is rarely efl^ective due to expressions swell. Much of the 
examples treated here and in [6], using our algorithm Chgt- 
Coords, can not be treated with this approach. 

It seems that the problem lies in the fact that we can not 
separate the computation of C^[f] from that of C/[/] which 
contains, very often, big expressions. 



An other disadvantage of the above method is that we 
have to restart computation from the very beginning if the 
target equation is changed. In the next section, we propose 
our algorithm ChgtCoords to compute the transformation 
alone and in terms of differential invariants. These invari- 
ants are provided by Cartan method for a generic / which 
means that we have not re-apply Cartan method if the target 
equation is changed and a big part of calculations is generic. 
Furthermore, the computation of ip in terms of differential 
invariants reduces significantly the size of the expressions. 

4. USING CARTAN'S METHOD 

In this paper, differential invariants are obtained using 
Cartan's equivalence method. We refer the reader to [51 1131 
1161 [7] for an expanded tutorial presentation and application 
to second order ODE. When applied Cartan's method fur- 
nishes a finite set of fundamental invariants and a certain 
number of invariant derivations generating the differential 
field of invariant functions. 

Example 7. Consider the EPB {J'^{C,C), "^3)- The PDE 
system ^ reads 



'dp- f{x,y,p)dx 

dy — pdx 
1 da: 





"2 




)■: 














dy — pdx 
da; 



S{a) 



With det(5(a)) 7^ 0. In accordance with Cartan, this system 
is lifted to the linear Pfaffian system 

S{a,) uj = 5(a) loj 

defined on the manifold of local coordinates (x, a, x, a) . After 
two normalizations and one prolongation, Cartan's method 
yields three fundamental invariants (p = y' and a — a^) 



h = 



fyp Dx fp: 

2a 



(23) 



and the invariant derivations 

a dy 2 a dp 2' ''^ da' 



(24) 



where = + p-§^ + f{x,y,p)-§^. 

When dim(<S/) = 0, the additional parameter a can be 
(post)normalized by fixing some invariant to some suitable 
value. In this manner one constructs invariants defined 
on M (not depending on the additional parameter). 

Theorem 2 (Olver jT6]). //dim(5j) = 0, then there 
exist exactly m functionally independent specialized invari- 
ants • ■ ■ ,/m[/]. 

Note that the invariants /i [/],•■ ■ ,Im [/] are functionally 
independent if and only if d/i [/] A • • • A dim [/] 7^ 0. Note 
also that if the function / is rational, then the specialized 
invariants /[/] : M C are algebraic functions. In the 
sequel, we use the notation h-j-.-k to denote the differential 
invariant ■ ■ ■ Xj (li). 



4.1 Computation of ip 

Suppose that the ©-groupoid <Sj is zerodimensional. Then, 
according to the theorem [S] there exists m functionally in- 
dependent invariants Ft '.~ Ik[f], 1 < k < rn. This implies 
that the algebraic (non differential) system 

{J'i(x) = Ji,...,F™(x) (25) 

is locally invertible and has a finite number of solutions 

x = F-'(Ji,...,J,„). (26) 

, Im on the source function / 



The specialization of Ji 
yields 



= F-\h[f\, ...,/„[/])■ 



(27) 



Let C denote the (non differential) characteristic set as- 
sociated to the system (|25|) w.r.t. the elimination rank- 
ing {x, } y {h, ...,/„}. Thus, C describes the 
inversion (|26p . The most simple situation happens when 
deg(C) = 1. In this case, the necessary form of the change 
of coordinates </? is the rational transformation defined by C. 

Example 8. Consider the EPB (J^(C,C),$3) and the 
target equation Ej introduced by G. Reid \17^ 



X X"^ 



The following invariants are functionally independent 



-20—, /l;31 

ax^ 



8^^, h = —rr + 8 



^i;3 



aa;-' 4a;2 
3 a; — 48 y + 16 pa; 
2i3 ■ 



We normalize the parameter a by setting /i,23 
characteristic set C is 



-20. The 



_3_ 

'32 



_3_ 

512' 



+ — /l;31^/l + 



4096 



Ii'ili/. 



256 



/l;31 + 



^ r3 

4096 ^I'^i' 



which gives the sought necessary form of (p. As a byproduct 
we deduce that the symmetry group TSf — {Id}. 

Let us return to the general situation, that is when deg(C) 
is strictly bigger than 1. We have two cases. First, deg(C) = 
deg(5j) and then tp is the algebraic transformation defined 
by C. Second, deg(C) > deg(<Sj). In this case, to obtain 
the transformation ip}, we have to look for m other function- 
ally independent invariants such that the new characteristic 
set C has degree equal to deg(<Sj). 

Example 9. Consider the EPB (J^(C, C), 'I>3) and the 
target equation £j (number 8 in ^91) 

-II -3 , 

y =y +^y 

which the corresponding symmetry group is 

rSf = {{x,y)^{x, \y)\X^ = l}. 

One can verify that Ii, 7i;i3 and Ii-ixi, when specialized 
on the considered equation, are functionally independent. In 



this case, the associated characteristic set C is 

3/i;33 - 2/1^) y 



P = 



y 



3(/l;3 + 1) 



-hx + -/l;33K 

+ 3 



— o^l;3^ + 7-^l;33A 
t o 4 



r/ie degree of this set is equal to 6 which is different from 
the degree of the symmetry groupoid. 

However, if instead of the above invariants we consider 
the invariants Ki := 7i;233/^i;3i, K2 := /i;234//i;3i and 
K3 :— /i;23i//i;3i, we obtain 

J = ~Kiy, 



6 



Ks + Ki 



This characteristic gives the necessary form of ip since it has 
degree two. 

4.1.1 Heuristic of degree reduction 

In practice, one has to search the invariants giving the 
required degree in the algebra of invariants. However, this 
is not an easy task since this algebra can be very large (al- 
though it is algorithmic). For this reason we provide an 
important heuristic which enables us to obtain the desired 
invariants. This heuristic is explained in the following ex- 
ample. 

Example 10. Consider the Emden Fowler equation {P)) 
and the T>-groupoid of transformations "1>3. We have already 
computed the corresponding symmetry groupoid. The spe- 
cialization of the invariants Ii, /i;i3 and /i;i33 gives three 
functionally independent functions. As explained above, we 
obtain the following characteristic set computed w.r.t. the 



ranking ■ 



y y y X y h y h-^-j 

9 
4 

/l;3/l 







p 


- 1 






f 


:!< 


x^ 






+8- 



>- A;33 
2 



h 

2 



r, ll!i . 3 Jl;33 

+2 



1 /l;3 - 



h 



(28) 



8/1.3' +6/l;33/l 
- 8Jl;3" + 6/l;33/l ■ 

Comparing with the V-groupoid of symmetries HlOjl we de- 
duce that, in contrary to y, the degree of x must be reduced 
to one. This can be done in the following manner. First, 
observe that the Lie defining equations of ^3, more exactly 
Xp = 0, implies that Xi{x) — where Xi — is the in- 
variant derivation \24^ . Now, differentiate the last equation 
of the characteristic set, which we write as x^ — Ax -\- B, 
w.r.t the derivation Xi. We find A-ix -\- B-i — 0. The co- 
efficient of x in this equation, which is invariant, could not 
vanish (since it is not identically zero when specializing on 

the Emden-Fowler equation). Thus, x — or explicitly 



With K 



Ii 



9h'-8Ii.,, 



(29) 



The necessary form of the change of coordinates ip is then 
given by Ij29[) and the two first equations of II28\) . 

The above reasoning can be summarized as follows 



Procedure ChgtCoords 

Input : £f and $ such that dim(<Sj) = 

Output : X — <^(x) the necessary form of the change of 

coordinates 

1- Find m functionally independent invariants 
(/i[/], /™[/]) defined on M. 

2- Compute a char, set C of the algebraic system (|26|) . 

3- If deg(C) = 1 then Return C. 

4- Compute Sf with Rosenfeld-Grobner. 

5- WHILE deg(C) ^ deg(<S/) DO 

Reduce the degree of C. 
END DO 

6- Return C. 



5. THE SOLVER 
5.1 Precalculation of 

5.1.1 The first step : the adapted V-groupoid 

Let $1, . . . , $7 C Jr(C^, C^) denote the ©-groupoids de- 
fined in the table 1 above. It is not difficult to see that 
$1 C $3 C $5 and $2 C "I'4 C $6 and finally $5, C $7. 

Let A{£f, <1>) := dim (aut(f /) n "l)) where aut(f/) is the 
contact symmetry ©-groupoid of the second order ODE £/. 
Let d, := A{£f, $i) for 1 < i < 7. 

Definition 10 (Signature). The signature of £f is 

sign(£:/) — ((di, d3, ds), (d2, d4, dg), d?) . 

Clearly, (di < ds < ds < d7) and (d2 < d4 < de < d7). 
Recall that the calculation of theses dimensions does not 
require solving differential equations. We shall say that the 
signature sign(i£/) matches the signature sign(iSj) if and only 
if d7 = d7 and (si = si or S2 ~ S2) where si and S2 stand 
for (di, ds, ds) and (d2, d4, ds) resp. Two second order ODE 
£f and £f are said to be strongly equivalent if 

3$ G {$1, • • • , $7}, G ip*£f = £f, d{£f, $) = 0. 



Lemma 2. If £f and £j are strongly equivalent then their 
signatures match. 

Definition 11 (Adapted D-groupoid). AT>-groupoia 
<I> is said to be adapted to the ODE £f ifd{£f, <I>) = and 
$ is maximal among "I>i, ■ ■ ■ , <1>7 satisfying this property. 





Transformations 


Equation number according 
to Kamkc's book 




X — X, y — r}{x, y), 


1, 2, 4, 7, 10, 21, 23, 24, 30, 
31, 32, 40, 42, 43, 45, 47, 50 


*3 


X = X + C, y = r}{x, y), 


11, 78, 79, 87, 90, 91, 92, 94, 
97, 98, 105, 106, 156, 172 


*5 


X = S,(x), y = ri{x, y), 


Null 




2 = 4(2:, j;), y = y, 


81, 89, 133, 134, 135, 237 


$4 


2 = i(x,v), y = y + C, 


11, 79, 87, 90, 92, 93, 94, 97, 
98, 99, 105, 106, 172, 178 


*6 


X = S,{x,y}, y = ri{y), 


80, 86, 156, 219, 233 


<l>7 


2 = 4(2:, j;), y = v(x-.y), 


3, 5, 6, 8, 9, 27, 44, 52, 85, 
95, 108, 142, 144, 145, 147, 
171, 211, 212, 238 



Table 1. 



The above table associates to each equation in the third col- 
umn its adapted groupoids. For instance, the first Painleve 
equation (number 3) appears in the last row which means 
that its adapted XJ-groupoid is the point transformations D- 
groupoid <1>7 . To the Emden-Fowler equation, number 11, 
we associate the O-groupoids $3 and $4. In the case of 
homogeneous linear second order ODE (e.g. Airy equation, 
Bessel equation, Gaufi hyper-geometric equation) we prove 
that, generically, the adapted D-groupoid is $4. 

5.1.2 The second step 

Once the list of adapted O-groupoids $ is known, we pro- 
ceed by computing the necessary form of the change of co- 
ordinates (/p G "1> using ChgtCoords. Doing so, we construct a 
Maple table indexed by Kamke's book equations and where 
entries corresponding to the index Sj are: 

1- the signature of f j, 

2- the list of the adapted O-groupoids <1> of £j, 

3- the necessary form of the change of coordinates iy3 G $. 
For instance, the entries associated to Rayleigh equation 
y + y + y — are: 

1- the signature ((0, 1, 1), (1, 1, 1), 1), 

2- the ©-groupoid <l>i, 

3- the necessary form of the change of coordinates 



P 

X 



-36- 



72 + 72 Ji + /2;i 
-1 



2^' 



559872 J|.i 



(4i + 2I6/1 J|;i + 216J|;i + 373248 



+15552/|.i + 11197447? -f 31104/i 
-1-1119744/1 +373248/? 



15552/? /|.] 



Procedure Newdsolve 
Input : £f 

Output : An equation £j in Kamke's book and the 
transformation tp such that ip,{£f) = £j 

1- Compute the signature of £/. 

2- Select from the table the list of equations £f such that 
sign(f j) matches sign(f/). 

3- FOR each equation £f in the selected list DO 

(i) Specialize, on £f, the necessary form of the change 
of coordinates associated to fj. We obtain ip. 

(ii) If ^ G $ and ipt{£f ) = £j then return (£j, ip). 
END DO. 



(30) 



with the normalization /2//2;i = 1. Invariants here are those 
generated by (|23p and H24I) plus the essential invariant x = x. 

5.2 Algorithmic scheme of the solver 

To integrate a differential equation £f our solver proceeds 
as follows 



It is worth noticing that the time required to perform steps 
(i)- (ii) is very small. In fact, it is about one hundredth 
of a second using Pentium(4) with 256 Mo. The second 
feature of our solver is, contrarily to the symmetry methods, 
neither the table construction nor the algorithm of the solver 
involves integration of differential equations. 
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